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Abstract The main difference between certain spectral problems for linear 
Schrodinger operators, e.g. the almost Mathieu equation, and three-term re- 
currence relations for orthogonal polynomials is that in the former the index 
ranges across Z and in the latter only across Z"*". We present a technique that, 
by a mixture of Dirichlet and Taylor expansions, translates the almost Math- 
ieu equation and its generalizations to three term recurrence relations. This 
opens up the possibility of exploiting the full power of the theory of orthogonal 
polynomials in the analysis of Schrodinger spectra. 

Aforementioned three-term recurrence relations share the property that their 
coefficients are almost periodic. We generalize a method of proof, due origi- 
nally to Jeff Geronimo and Walter van Assche, to investigate essential support 
of the Borel measure of associated orthogonal polynomials, thereby deriving 
information on the underlying absolutely continuous spectra of Schrodinger 
operators. 
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1 The almost Mathieu equation and orthogonal polynomials 

The point of departure of our analysis is the almost Mathieu equation (also known as the 
Harper equation). We seek A G M and {an}n(^z S that satisfy 

On-i — 2Kcos(an + /3)a„ + a„+i = Aa^, n £ Z. (1-1) 

a, f3 and n ^ being given real constants. 

The almost Mathieu equation features in a number of applications Q and has been 
already extensively studied |||, ^ |l^, 15 1. The purpose of our analysis is not to reveal 
new features of the spectrum of ( |1.1| ) per se, since the latter is quite comprehensively 
known. Instead, we intend to demonstrate that the almost Mathieu equation exhibits an 
intriguing connection with orthogonal polynomials, a connection that lends itself to far 
reaching generalizations. 

Let a; = e'", bi = Ke'^^ and 62 = «^e~'^ = 61. We rewrite ( |1.1D as 

(6iu;2'' + 62)an = '^"(an-i -Aa„ + o„+i), neZ (1.2) 
and consider the Dirichlet expansion 



The choice of a specific branch of the square root of 61 is arbitrary. Note that, since 
\uj\ = 1, it is easy to demonstrate that {an}n£Z G implies convergence of the series 
for all z e C 0, 0]. 

We multiply ( |l.2D by exp |6^u;"zj' and sum for n G Z. Since 

J 00 1 00 1 

y'iz)=b! «na;"e''?-"^ y"(z) = 61 ^ aW^"^ 



we readily deduce that y obeys the functional differential equation 

y"iz) + b2yiz) = bf^ {u;-Wi^-'z) - \y'{z) + uy'icoz)] . (1.3) 



The solution of ( |1.3| ) is determined uniquely by the values of y{0) and y'{0). 

Dirichlet expansions have been employed by Gregori Defel and Stanislav Molchanov 
Q to investigate the simplified spectral problem 

a„_i - 2Kg"a„ + a„_i = Aa„, n G Z, 

where k G ]R\{0} and q > 0, and they have demonstrated that it obeys another functional 
differential equation. Both the Derfel-Molchanov equation and (|1.3D are a generalization 



of the pantograph equation, that has been extensively analysed in [1^ and |1C]. However, 
it is important to emphasize that ( p,.3D has an important feature that sets it apart from 
other functional equations of the pantograph type, namely that, unless u; G M, its evolution 
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makes sense only for complex z and proceeds along circles of constant \z\, emanating from 
the origin. 

Inasmuch as the equation (|1.3| ) can be analysed directly, our next step entails expanding 
it in Taylor series. Thus, letting 



m=0 



' ^ 'Tin I 

substitution in (|l.3D readily yields 

_ 1 

ym+i = Vb[ ^ (2 cos am - A)ym - 62ym,-i, m = l,2, (1.4) 

It is beneficial to treat Um as a function of the spectral parameter A and to define 

ym{t) ■■= h^'^'\m{-{hih2f'h), m = 0, 1, . . . . 
Brief manipulation affirms that (|l.4| ) is equivalent to 

ym+i{t) = {t+'^cosaTn^yrn{t) -yrn-i{t)-, m = 0, 1, (1.5) 



To specify the solution of (1.5) in a unique fashion we need to choose yo and yi, which, of 
course, corresponds to equipping ( |1.3| ) with requisite initial conditions. We note in passing 
that t = —kX, hence real initial values in ( |1.3D correspond to yo,yi S 

Each solution of ( fl.Sj) is a linear combination of two linearly independent solutions. 
Setting a = 2/ K, we let 

r_i(t) = 0, 

ro{t) = 1, 

rm+i{t) = {t + a COS am)rm{t) - rm-i{t), m = 0, 1, . . . , (1.6) 

and 



s_i(i) = 
so{t) = 1 

Sm+i{t) = {t + a COS a{m + l))sm{t) - Sm-i{t), m = 0, 1, (1.7) 

It is trivial to verify that {rm}m&Z+ and {sm+i}meZ+ are linearly independent, hence 
they span all solutions of ( |1.5| ). Moreover, we observe that each rm{t) and Sm{t) is an 
mth degree monic polynomial in t. This is a crucial observation, by virtue of the Favard 
theorem [^: Given any three-term recurrence relation of the form 

P-i{t) = 0, 
Vo{t) = 1, 

Pm+l{t) = it + Cm)Pm{t) - dmPm-l{t), m = 0, 1, . . . , 



4 



A. Iserles 



the monic polynomial sequence {pm}meZ+ is orthogonal with respect to some Borel mea- 
sure dip, i.e. 

Pm{t)Pn{t) d(p{t) = 0, m^n. 



The essential support H of d^p is of central importance to the matter in hand, since it is 
an easy consequence of, for example, the nth root asymptotics of orthogonal polynomials 



[ 16| that, as long as H is bounded, the sum Ylm=QPm{t)z'^ /m\ converges for z G C, 
whereas convergence fails for z 7^ for i S. Thus, travelling all the way back from 
orthogonal polynomials to functional equations, and hence to the almost Mathieu equation, 
we deduce that, subject to the linear transformation A = —tjn, essential supports of the 
Borel measures corresponding to ( p..6D and ( |1.7D result in the essential spectrum of ( |1.1| ). 

We note that, in principle, the Favard theorem falls short of producing a unique mea- 
sure and it is entirely possible that there exist many Borel measures that produce an 
identical set of monic orthogonal polynomials. This, however, is ruled out by the deter- 
minacy of the underlying Hamburger moment problem and the latter can be affirmed for 
both ( |1.6| ) and ( [1.7] ) by the Carleman criterion [0. 

The remainder of this paper is devoted toward the determination of H. In Section ^ 
we demonstrate, by extending a technique due to Jeff Geronimo and Walter van Assche, 
that the essential support - and, indeed, the underlying Borel measure - can be specified 
explicitly when a/vr is rational. We prove in essence that d(/9 is a linear combination of 
Chebyshev measures of the second kind, supported on a set of disjoint intervals. 

Various results on irrational a/vr are reported in Section |3|. Inasmuch as the general 
form of d(/9 is currently a matter for conjecture, we derive a number of results that, beside 
being of interest on their own merit, are fully consistent with known results about the 
almost Mathieu operator, in particular with the theorem that the absolutely continuous 
spectrum of (|l.l| ) is (for irrational a/vr) a Cantor set. 



The equation (1.1) has been extensively studied in the past and, inasmuch as there 
are few outstanding conjectures, our knowledge of the spectrum of the almost Mathieu 
operator is quite comprehensive. Although the approach of this paper introduces a new 
perspective, there is no claim that Sections add to the current state of knowledge 
of Schrodinger spectra. This state of affairs is remedied in Section ^, where the frame- 
work of our discussion undergoes a far reaching generalization. Firstly, we demonstrate 
that general periodic potentials with a finite number of Fourier harmonics lend them- 
selves to similar analysis, except that, instead of orthogonal polynomials, the outcome is a 
generalized eigenvalue problem for a certain matrix pencil. Secondly, we prove that a mul- 
tivariate extension of the almost Mathieu equation can be 'transformed' by our techniques 
to a problem in (univariate) orthogonal polynomials. 

A possible application of our analysis, and in particular of Section ^, is numerical 
computation of the essential spectrum of ( |1.1| ) and of its generalizations. We do not 
pursue this further in the present paper. 

The observation that ( |l.lD is 'almost' a three-term recurrence relation - only the in- 
dex range is wrong - hence that the almost Mathieu equation might be connected with 



orthogonal polynomials, is not new. Most prominently, it had been made in |12], where it 
motivated a very interesting generalization of Chebyshev polynomials. The most impor- 
tant innovation in the present paper is in a technique that reduces the index range from Z 
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to Z"*" by passing from Dirichlet to Taylor expansions and which can be extended to cater 
for a substantially more general problem. 

2 Orthogonal polynomials with periodic recurrence coeffi- 
cients 

The focus of our attention in this section is the three-term recurrence 

P-i{t) = 0, 
Po{t) = 1, 

Pm+l{t) = {t - am)Pm{t) - Pm-lit), m = 0, 1, . . . , (2.1) 

where the sequence {am}meZ+ is i^-periodic, 

am+K = am, m = 0, 1, (2.2) 

Note that both (|]|) and assume this form when a/vr is rational. Our objective is 
to determine the Borel measure that renders {Pm}m&Z+ i^^o an orthogonal polynomial 
system (OPS). 

In Q the authors consider the following problem. Let {Qm}m.e'Z+ be an OPS whose 
measure has an essential support Hq ^ [—1,1] and let T be a given A^th degree polynomial. 
Setting H = T~^{Eq) (the latter set is, generically, a union of < disjoint intervals), 
they derive a new OPS, whose Borel measure is supported in S, explicitly in terms of 
{Qm}meZ+ ■ This construction is intimately related to the discussion of this section, except 
that we need, in a manner of speech, to travel in the opposite direction. As it turns out. 



the Borel measure associated with (2.1) inhabits a sets of disjoint intervals and we identify 
it by choosing an appropriate polynomial transformation T. 
Let 

Qn{t) ■■=P(n+l)K-lit), n = 0,l,.... 

Note that g_i = and, moreover, (p.l|) and (|2.2[) imply 



PnK{t) = {t - ao)qn-i - PnK~2it), n=l,2, .... (2.3) 
We seek polynomials ai, (3i, i = 0,1, . . . , K — 1, such that 

PnK+l>{t) = ai{t)qn-l{t) -bi{t)pnK-2{t), ^ = 0,1, . . . , K - I, 71=1,2,.... 

Because of ( |2.3D and the definition of g„, we have 

a^i{t) = 1, 6_i(t) = 0, 

ao(t) = t-ao, h{t) = 1. ^ ' ' 

We next substitute in the recurrence relation ( p.lj ) and, by virtue of (p. 21), obtain 

PnK+l+l{t) = {t - ae+i)pnK+e{t) - PnK+e-l{t) 

= (t- ae+i){ai{t)qn^i{t) - bi,{t)pnK-2{t)} 
- {ai_i{t)qn-i{t) - be-i{t)pnK-2{t)}. 
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Thus, comparing coefficients, we derive tlie recurrences 

ae+iit) = {t - ai+i)ai{t) - ai^i{t), 

be+i{t) = {t-ae+iMt)-be_i{t), £ = 0,1, 

wliich, in tandem witli ( p. 4) ), determine {(ie,bi}fjQ^ . 
Let i = K - 1, tlien 



.K-2, 



(2.5) 
(2.6) 



PnK-2{t) = 

and, sliifting the index, 

P{n+l)K~2it) 

Substituting both expressions into 



aK-iit)qn-i{t) - qn{t) 
bx-iit) 

_ aK-i{t)qn{t) - qn+ijt) 
bx-iit) 



P{n+l)K-2{t) = aK~2{t)qn~l{t) - bx ~2{t)PnK ~2{t) 

yields the recurrence relation 

qn+lit) = {aK-l{t) - bK-2{t))qn{t) " A„_2 (O^n-l (*) ) 

where 



(2.7) 



Af(t) = det 



be{t) be+i{t) 



0,1,...,K-1. 



We multiply (|2.6| ) by ae{t), (|2.5| ) by be{t) and subtract from each other. This readily 
affirms by induction that 

A<.(t) = A^_i(t) = ••• = ! 

and ( p.?! ) simplifies into 



^+i(t) = (ax^iit) - bK-2{t))qn{t) - q-n^iit), n = 0, 1, 



(2.8) 



Note that (|2.8| ) is consistent with n = 0, since q-i = 0. To further simplify the recurrence, 
we observe that qo{t) = aK-i{t) — bK-2{t), hence, letting 



qn{x) :-- 



.it) 



qoit)' 



n 



-1,0 



, w, . . . , 



where x = qoit), we obtain the three-term recurrence 

g_i(x) = 0, 
qo{x) = 1, 
qn+i{x) = xqn{x) - qn-i{x), n = 0, 1, 

Thus, each ijn is an nth degree monic polynomial and, by virtue of the Favard theorem, 
{qn}n£'L+ is an OPS. It can be easily identified as a shifted and scaled Chebyshev polynomial 
of the second kind, 

g„(x) = 2"C/„ (ix) , n = 0,l,.... 
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We thus deduce that 



qn{t) = 2''qo{t)Un {^oit)) , n = 0, 1, . . . . 



(2.9) 



Before we identify the underlying Borel measure, let us 'fill in' the remaining values of 
Pm- By definition, PnK-i = Qn-i, Pn{K+i)-i = Qn, hence the recurrence (2J) gives 



(x - ai )pnK it) - PnK+l (t) = 

-PnK+e-lit) + (t- ai+i)pnK+eit) - PnK+e+lit) = 

-P{n+l)K-3{t) + it- aK~l)p(n+l)K-2{t) = 

This is a linear system of equations, which we write as 

Ax-lPn = fin, 

where 



qn~i{t), 

0, i 

qn{t)- 



1,2,. 



(2.10) 



t — ai —1 
— 1 t — a2 



-1 



-1 t — as 



-1 t-am-i -1 

— 1 t — Or 



m = l,2,...,K-l., 



Pn 



PnK{t) 
PnK+l{t) 

P{n+l)K-3{'t) 
P(n+l)K-2i't) 



and 



qn 



qn-i{t) 







2n{t) 



We expand the determinant of Am in its bottom row and rightmost column. This results in 
a three-term recurrence relation and comparison with (2.4) and (2.6) affirms that det Am = 
hm{t). Hence, solving ( 2.1CD with Cramer's rule, we deduce that there exist (ET — 2)-degree 
polynomials and 6^, ^ = 0, 1, . . . , — 2, such that 



PnK+l{t) 



ai{t)qn-i{t) + bi{t)qn{t) 



^ = 0,1, 



hK-i{t) 

Bearing in mind the definition of ai and bi, 

PnK+£(t) = ai{t)qn-l{t) - bg{t)pnK -2{t) , 

we obtain from ( 2. Ill ) the identity 

ai){t)qn-i{t) + bi,{t)qn{t) bi{t){aK-2{t)qn~2{t) + bK-2{t)qn-i{t)) 

- = a£[t)qn-i[t) 



(2.11) 



bK~i{t) 



bK-iit) 
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We next substitute 

qn-2{t) = (ax-lit) - hK-2{t))qn-l{t) - qn{t) 

[pace (|2.8| )) and rearrange terms, whereby 



{h{t) -h{tyaK-2{t)}qn{t) 
= {-CLi + ae{t)bK-i{t) - {aK-i{t)bi{t) - hi{t)hK-2{t))aK-2{t) - bi{t)bK-2{t)}qn-i{t). 

However, consecutive orthogonal polynomials qn-i and g„ cannot share zeros [||, therefore 
both sides of the last equality identically vanish and we derive the explicit expressions 

de{t) = ae{t)bK-i{t) - {aK-i{tMt) - be{t)bK-2{t))aK-2it) - be{t)bK-2it), (2-12) 
kit) = be{t)dK-2{t). (2.13) 



Letting £ = K -2m gives bK~2{t) = bK-2(t)dK-2{t) and we substitute this into 

( 2.12 ). The outcome is 



In particular. 



de{t) = ae{t)bK-i{t) - aK-i{t)bi{t)aK~2{t)- 
K — 2 and the definition of result in 



(2.14) 



(1 + aK~l{t)bK~2{t))dK-2{t) 



aK-2{t)bK-l{t) = aK-l{t)bK-2{t) + ^K~2{t) 
1 + aK-i{t)bK-2{t). 



Since aK-ibK-2 ^ — 1, we conclude that 0^-2 = 1 and substitution in ( 2.13 ) and ( p.l4| ) 
yields the explicit formulae 



det 



ai{t) aK-i{t) 
h{t) bK-i{t) 



k{t) = be{t), e = 0,l,---,K -2. 



Theorem 1 The OPS {pm}meZ+ has an explicit representation in the form 

1 



PnK+e{t) 



bK-i{t) 



det 



ai{t) aK~i{t) 
be{t) bK-i{t) 



In-l 



(t) + bi{t)qn{t) 



(2.15) 



where n = 0, 1, . . ., i = 0,1, . . . , K — 1 and the OPS {qn}n(^Z+ satisfies the three-term 
recurrence (^^). □ 



Note that letting £ = K — loic£ = — lin ( p. 15 ), in tandem with ( |2.4| ), results in 
P{n+i)K-i = qn and PnK-i = respectively, as required. 

Let t G H. Then, by the discussion preceding the representation ( |2.9D , we know that 
^qo{t) G [—1,1], and there exists 6 G [— tt, vr] such that ^qo{t) = cos 6*. Since, by the 
definition of Chebyshev polynomials of the second kind. 



C/„(cos 9) 



sin(n + 1) 



smf 
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( |2.9| ) implies that 

Substitution into ( |2.15| ) results in 
2No(i) 



sin(n + 1)0 



sm ( 



PnK+e{t) 



bx-iit) sin 6 



det 



bi{t) bK-l{t) 



sinnO + ^^(t) sin(n + 1)^ 



We next proceed to determine the essential support H of the Borel measure dip which 
corresponds to the OPS {pn}nGZ+- we have already mentioned, this is very similar to 
the construction of Geronimo and van Assche in with {qn}neZ+-: ^Qo and [—1, 1] plying 
the role of {Qn}neZ+ ^ ^ and Hq respectively. This similarity notwithstanding, there are 
some important differences - not least that our argument advances in an opposite direction 
to that of Geronimo and van Assche - and we present here a complete derivation of H. We 
commence by observing that, by virtue of Theorem 1, everything depends on the support 
of qn{x{t)) = 2^Un {^qo{t)^ , n = 0, 1, . . .. Thus, we seek a Borel measure dip such that 



In,; 



qn{x{t))qm{x{t)) dil){t) =0, n,m = 0,1, 



n ^ m. 



It follows from the defintion of g„ that 



-tn+k 



Un {\qo{t)) Um {\qo{t)) dm- 



Similarly to |^], we seek the inverse function to x = \qo{t). Let ^^i < < ■ ■ ■ < be all the 
minima and maxima of go iii 1^ (of course, s < K — 2) and define ^0 = —00, ^s+i = 00. In 
each interval [.^j, .^j+i], j = 0, 1, . . . , s, the function ^qo{t) is monotone, hence it possesses 
there a well-defined inverse. We denote it by Xj{x), hence ^qo{Xj{x)) = x. Changing the 
integration variable, we have 



^n+mj- 1^'^' Un {\q^{t)) Um {hoit)) dm 



j=0 



Unix)Um{x)dlPiXjix)). 



(2.16) 



We recall that {Un}neZ+ is an OPS with respect to the Borel measure (1 — x^)2 dx, 
supported by x G [—1, 1]. Thus, for every j = 0, 1, . . . , s we distinguish among the following 
cases: 

Case 1: go(Ci) < -2 and 2 < go(Cj+i)- 

We stipulate that dtjj{Xj{x)) vanishes for all 



X G 



ko(e,),ko(e,+i)l 



Since Xj increases monotonically in the contribution of this interval to ( p.l6| ) is 

\n{x)Um{x)d^{\X,{x)\). (2.17) 



10 



A. Iserles 



Case 2: qo{Cj+i) < -2 and 2 < qo{Cj)- 

Likewise, we require that the support of dip{Xj) is restricted to [—1, 1]. qo decreases mono- 



tonically within and straightforward manipulation affirms that (2.17) represents 

the contribution of this interval to ( |2.16| ) . 

Case 3: min{go(Cj), 9ofe+i)} > -2 or max{go(Cj), go(Ci+i)} < 2. 

In that case we cannot fit [—1,1] into hence we stipulate that d'ip{Xj) is not 

supported in 

Let z^i < z^2 < ■ ■ ■ < i^r be all the indices in {0, 1, . . . , s} such that either Case 1 or 
Case 2 holds. We require that r > 1. Then (|2.16D reduces to 



Un{x)Um{x)Y. dlPi\X,,ix)\). 

Since the Hamburger moment problem for the Chebyshev measure of the second kind is 
determinate, it follows that necessarily 

r 

dV'(|X,,(x)|) = (1 - x2)i dx, X G [-1, 1]. 



Theorem 2 The orthogonality measure corresponding to the OPS {Pm}meZ+ is supported 
by 

where for each i = 1,2, ... C [S^u^^^u^^^] is the unique interval such that |(7o(OI = 2 
at its endpoints. 

Proof. Follows at once from our construction. □ 

Figure 1 displays two examples of the present construction, for different cases of q^. In 
each case H is the union of the 'thick' intervals. 



Harking back to (|1.6| ) and (|1.7|) , we let am = —a cos am and am = —acos{m + l)a, 
m = 0,1,..., respectively, where a = IttL/K. Thus, {am}meZ+ is indeed i^-periodic. 
Unsurprisingly, the outcome of our analysis are the familiar spectral bounds ^. The 
merits of our approach are, however, not just in providing an alternative proof of known 
results but also in extending the framework to the multivariate case in Section ^ 

We mention in passing that the analysis of this section can be easily extended to 
recurrences of the form 

P^iit) = 0, 
Po{t) = 1, 

Pm+lit) = (t - am)Pmit) - (3mPm-l(t), m = 0, 1, . . . , 

where both {am}m&Z+ ^iid {(3m}mS:I,+ i^-periodic. This, however, is of little relevance 
to the theme of this paper. 




Figure 1: The sets E for two different polynomials go- 
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3 Orthogonal polynomials with almost periodic recurrence 
coefficients 

The three-term recurrence relations ( |1.6D and (L7) assume, a/vr being irrational, almost- 
periodic recurrence coefficients and this state of affairs is even more important in a mul- 
tivariate generalization of (O) in Section ^. Unfortunately, no general theory exists to 
cater for orthogonal polynomials with almost periodic recurrence coefficients. The theme 
of the present section is a preliminary and - in the nature of things - incomplete investi- 
gation of the case when an irrational a/vr is approximated by rationals. In other words, 
we commence with the X-periodic recurrence ( |2.lD , except that we will allow the period 
K to become unbounded. 

The motivation for our analysis is an observation which is interesting on its own merit. 
Denote by ax the value of qx-iiO) for ae = — cos i = 0,1, . . . , K — 1, i.e. 



ax = det 



cosf 1 
1 cos 



4-77 
K 



COS 



2{K-2)-K 
K 



COS 



1 

2{K-1)tt 
K 



K = l,2,.... 



Computation indicates that 



<T4L = 0, 



lim CT4L+2 
L— >oo 



hm <T4L+i 

L^oo 



lim (T4L+3 

L^oo 



2^/3 



It is easy to prove that a^L = for all L > 1. Thus, let K = AL and set £ = L in 



Vi = cos - ^^£-2, 

K 



1,2,...,4L-1 



(3.1) 



(3.2) 



(note that (|3.2[) , in tandem with = 0, vq = 1, yields ax = vk-i)- This yields 
+ 1^1-2 = and we claim that, in general. 



t'L+k + (-l)*''j^L-fc = 0, 



k 



-1,0,...,L-1. 



(3.3) 



We have already proved ( |3.3| ) for k = and it is trivially true for k = —1. We continue 
by induction and assume that ( p.3[ ) is true for k = — l,0,...,s — 1. Letting £ = L it s in 
(|3.2| ), we have 



^L+s = — sm 27;^'L+s-l - 

l^L-s = +Sm.^UL^s~l- 'I^L-s-2- 

We multiply the second equation by (— l)'^ and add to the first, thus 

{VL+S + {-lyi^L-s) = -Sinfg(l/L+,_1 + {-ly-^VL-s+l) - {vL+s-2 + {-1)'''^ VL-s+2) 

and ( ^.31 ) follows at once. 
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Hence, letting A; = L — 1 in ( |3.3| ) and recalling that z/_i = 0, we obtain z/2l_i = 0. 
Moreover, similarly to (|3.3| ), we can prove that 

l^SL+k + (-l)V3L_fc = 0, /c = -1, 0, . . . , L - 1 

and k = L — 1 gives 

0-4L = i^4L-l = {-l)^V2L~l = 0. 

This completes the proof of the lemma. □ 



Other observations in ( |3.1| ) are also true and in the sequel we prove them in a gener- 
alized setting. 

Given a G (—1, 1), we define 



A-i{t) = 0, AQ{t) = l, 

An{t) = tinAn-l{t) - An-2{t), n = 1, 2, . . . , 



(3.4) 



where {in]^=i is a given real sequence. To emphasize the dependence on parameters, we 
write, as and when necessary, A„( • ) = An[ • ; ^i, • • • > in)- 
An alternative representation of An is 



An{t) = det 



til 1 

1 t6 



■•. 1 

1 tin 



n = 1,2,..., 



hence An is an nth degree polynomial. We observe that 

A„(0) = 



{-ly 





: n = 2s, 
: n = 2s + 1. 



Moreover, differentiating with respect to t, we obtain 

til 1 



<(t) = 5]det 
fc=i 



1 



ik 

1 tik+i 



1 t^r, 



therefore, expanding in the kth row, we derive the identity 



^n(*; • • • 1 '^n) — X! Cli • • • 1 '^fe-l)^n-fc(*; Cfc+li 

fe=l 



(3.5) 
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Proposition 3 An"^ (0) = whenever n + r is odd, hence the polynomial A^ has the same 
parity as n. 

Proof. By induction on r. The assertion is true for r = 0. Moreover, repeatedly 
differentiating ( ^.51) with the Leibnitz rule and letting t = 0, we obtain 

2r / \ In 
£=0 \ * / fc=l 



But 



{k — 1) -\- 1 is even 



(2n -k) + (2r - ^) is odd. 



therefore for all ^ = 0, 1, . . . , 2r and k = 1,2, . . . ,2n the induction hypothesis affirms that 
at least one of the terms in the product vanishes. Similar argument demonstrates that 

4nil(0)=0. □ 

We therefore let for all n = 0, 1, . . ., s = 1, 2, . . ., 



n+s— 1 



r=0 
n 

{2r+l),2r+l 



^2n+l(i; 6> • • • ) 6n+s) — E -^2n+l,s* 

r=0 

Substitution into (|3.4D (where we replace by ^n+s_i) results in the recurrences 



B. 



(2r) 
2n,s 

(2r+l) 
2n+l,s 



^ n(2r-l) „{2r) 

?s+2n-l-D2n-l,s ^2n,-2,s' 

?s+2n-D2n,s -"2n-l,s- 



(3.6) 
(3.7) 



Given g € C, we recall that the q-factorial symbol is defined as 

n-l 

{z;q)n=l[il-q''z), z£C, nGZU{oo}, 
whereas the q-binomial reads 



n 
m 



< m < n. 



Lemma 4 Let = + a'', s = 1^ 2, . . ., where g G C is given. Then, for every 
n = 0, 1, . . ., s = 1, 2, . . . and r = 0, 1, . . . , n, 



S^''^ = (-l)"+'-g-K2"+«-^-^) 

n(2r+l) _ / ^.n+r-(r+^){2n+s~r) 



n + r 
2r 



(3.8) 



n + r + 1 
2r + 1 



-'Z"+^-'^;g)2r+i. (3.9) 
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Proof. By induction on n, using (3.8) and ( |3.9| ). Obviously, the assertion of the lemma 
is true for n = 0. Otherwise, for even values, 



q ^ -\- q ^1 02n-l -D2ri-2,s 



( 

l^_l^n+r ^-{r-^){2n+s-r-l) ^ 



(-1) 



n+r—l—r(2n+s—r- 



qn+ 2 _|_ g " 2 



n + r — 1 
2r 



n + r — 1 
2r - 1 



2r-l 



n+s—r—1. 



-,q)2r 



I ^\n+r „~r{2n+s~-r~};) (^i g)n+r-l ( n+s-r.\ 
l-J-J ^ ^ ^'l N 7 ^ \-1 ,q)2r-l 

X {(1 + g2n+^-l)(l - q^') + _ + 

(_l)"+''-^-'-(2n+«-r--|) 



n + r 
2r 



n+s— r . 



This accomplishes a single inductive step for (3.8). We prove ( p. 91 ) in an identical manner, 
by considering odd values of n. □ 

Recall that our interest in the polynomials An has been sparked by the observation 
(|3.1| ). Thus, we require to recover cosine terms, and to this end we choose q of unit 
modulus. 

Proposition 5 Suppose that g"^ = 1 and 7^ 1 for m = 1,2, ... ,n. Then -B2n+i 1 ~ 
foranO<r<^. 



Proof. Let 2r < n. We have from (|3.9|) that 



B. 



(2r 
2n+l 



+1) _ ( -,\n+rjr+i:)ir+l) (g; g)n+r+l f 



However, 



(9;9)n 



+r+l 



{q; q)n-r 

whereas, because of our restriction on r. 



{q;q)2r+i{q]q)n-r 



n (i-g')=o, 



2r+l- 



2r+l 

{.q\q)2r+i= n(l-«')/0' 
^=1 



since g is a root of unity of minimal degree n + 1. The proposition follows. □ 

2mm 
n+1 



Corollary Let q = exp ^^xr, where m and n are relatively prime. Then it is true that 



lim A2n+i{t) = 

n^oo 

for every t € ( — 1, 1). 

Proof. Straightforward, since A2n+i{t) = O^^t^"^. □ 



(3.10) 
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Proposition 6 Suppose that uj = q2 is a root of unity of minimal degree 2n + 1. Then, 



for all r = 0, 1, . . . , n it is true that 



where (p = argw. 

Proof. Since = 1, it follows from (U) that 



{-'-) ^2n,l - 1 ^ ^'-( ^ 7 ^ \~1 ^iQhr- 

[q;q)2r[q;q)n-r 

But 

iq;q)n+r{q~''^^;q)2r ^ yr _ ^2^+1^ 
(q; q)n-r ^tZr 

Moreover, [q] q)2r 7^ for r = 0, 1, . . . , n, since 2n + 1 is the least nontrivial degree of the 
root of unity q, and we deduce that 

But 

1=1 1=1 1=1 

and, likewise, 

n (1 - q^'"-^) = n '^''+'(^''+' - ^"''"') = (-1)'^ n('^''+' - ^-''-')'. 

l=-r i=-r e=0 



Consequently, 



n R(2r) _ nLo(^ 



r-l/,,,2£+l _ ,.,-2e-l\2 



This is precisely the identity ( |3.11| ). □ 

Next, we consider progression to a limit as n — > 00 - ( ^.ip is a special case. Thus, 
suppose that we have a sequence $ = {(j)n}n€i, where 0„ = 27rmn/{n + 1), m„ G Z+, 
X C Z"*" is a set of infinite cardinality and 

lim 4>n = (j) € [0,27r). 

n— >oo 

Set 

C{t, = Jim (-l)["/2U„(t;ei"\e?\ . . . 
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where ^^^^ = 2cosicl)n, ^ = 1,2, ...,n, n = 1,2,.... Thus, ^^"^ = qn + Qn^, where 
qn = exp Consequently, according to Proposition 5, if X consists of only odd indices. 



necessarily C(t, $) = 0. This proves, icidentally, that cj2l ^ in (3.1 



2 

Lemma 7 Suppose that Z C and that = 0. Then, provided that m„ = o(n3), it is 
true that C(t, = (1 - t^)"! . 



Proof. Since 



sin 2£(/)„ 21 
we deduce from ( p.ll| ) that 

^2r 



+ 



r 

According to the assumption, 0[n(j)'^) = o(l) and the lemma follows from 

^ /2r\ t^ _ 1 

□ 



Letting t = h affirms the remaining part of (3.1) 



Similarly to the last proposition, it is possible to derive an explicit expression for 
C{t,^), provided that cp/ir is rational and that n((/> — (pn)^ = o{l) as n ^ oo. The 
derivation is long and it will be published elsewhere. It suffices to mention here the 
remarkable sensitivity of C{t, <I>) to both the choice of X and to the specific nature of (f). 

What happens when (p/n is irrational? This is, as things stand, an open problem. It 
is possible to show that, formally. 



i=0 k=l 2/c<A {q;q)i' 

where q = e^''^. The latter series can be summed up by means of the Gaui5~Heine theorem 
Q for \q\ < 1 and, after simple manipulation, for \q\ > 1. Unfortunately, because of a 
breakdown in Holder-continuity across \q\ = 1, it is impossible to deduce its value on the 
unit circle from the values within and without by means, for example, of the Sokhotsky 
formula Q. 

Clearly, there is much to be done to understand better the behaviour of the p^s when 
the period K becomes infinite. In this section we have established few results with regard 
to the values at the origin. They should be regarded as a preliminary foray into an 
interesting problem in orthogonal polynomial theory cum linear algebra and we hope to 
return to this theme in the future. 



4 Generalizations 

There are two natural ways of generalizing an almost Mathieu equation ( |l.lD , by either 
specifying a more general periodic potential or replacing the index by a multi-index. Re- 
markably, the basic framework of this paper - replacing a doubly-infinite recurrence by a 
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functional equation which, in turn, is replaced by a singly-infinite recurrence - survives 
both generalizations! In the present section we describe briefly this state of affairs. 

Firstly, suppose that the cosine term in ( |1 . ![ ) is replaced by a more general harmonic 
term and we consider the spectral problem 

r m ^ 

ttn-i ~ 2 < ^ Ke cos{n£6 + ipi) ^ an + a„+i = Aa„, n € Z. (4-1) 



We assume that ki, ■ ■ ■ , i^m £ ^ and, without loss of generality, that ^ 0. Letting 
K| = e''^^K£, K*_(, = e-'^^^Ke, i = l,2,...,m, 



in 

q = e , we set 



and Kq = 0. Therefore (4.1) assumes the form 



■ 2m 



an-1 - J2 ^i-mQ""' «n + g'""a„+i = A(?™"an, n G Z. (4.2) 
U=o J 

Let c e C \ {0} and consider the Dirichlet series y{t) = J2'i^=-oo'^new{'^Q^i}- Since, 
formally, 

oo 

yW(t)=c^ Y: g-Ve^""*, £ = 0,1,..., 



we obtain from (4.2) the functional differential equation 

2m 

Y.4-mC"'''y^'\t) = q'^y^'^Hqt) - \y^^\t) + q~^y^'^\q-H). (4.3) 
i=0 



The derivation is identical to that of (O) and is left to the reader. 

In line with Section |l|, we next expand the solution of ( [4.3| ) in Taylor series, y{t) = 
J2'^=oPnt"' /nl. This readily yields 

2m 

J2 4-mC^''Pn+i = (^7"+™ + g-"-™ - X)Pn+m, n = 0, 1, . . . , 

1=0 

hence, replacing n + m by n, 

m 

K}c~^Pn+i = (2 COS n9 — X)pn, i = m,m + 1, . . . . 

£=—m 

Finally, we choose c = expjiV'm/'^}, hence Kj!j^c~™ = = Km € M\{0}. We thus define 
ai = c~^K^/Km, < m and replace A by —Xk^- This results in the recurrence 

m 

aipn+l = (A - Pn)Pn-, n = m,m + l,..., (4.4) 

e=-m 

where 

cos n9 

fin = , n = m,m + l, 



Schrddinger spectra and orthogonal polynomials 



19 



Note that, inasmuch as the a^s may be complex, we have = ai, £ = l,2,...,m, 
ao = 0. 

The recurrence (|4.4| ) is spanned by 2m hnearly independent solutions. However, unless 
m = 1, it is no longer true that, for appropriate choice of po,pi, . . . ,p2m-i, each pn is a 
polynomial of degree n + k for some independent of n. Indeed, it is easy to verify that 
the degree of pn increases roughly as [n/m]. Hence, orthogonality is lost. Fortunately, an 
important feature of orthogonal polynomials, namely that their zeros are eigenvalues of a 
truncated Jacobi matrix Q, can be generalized to the present framework. It is possible 
to show that the zeros of pn are generalized eigenvalues of a specific pencil of 'truncated' 
matrices and this provides a handle on their location. We expect to address ourselves to 
this issue in a future publication. 

Another generalization of (1.1) allows the index n to be replaced by a multi-index 



n = (ni, n2, . . . , n^) G Z"'. Thus, let G be the ^th unit vector, £ = 1, 2, . . . , d, and 
consider the spectral problem 



d / ^ \ 

^(On+ef + fln-ej " 2kCOS ^ a£n£ + /3 j 
£=1 \l=l 1 



«n = Aan, n e TJ^. (4.5) 



In line with Section |T|, we let 

6i=6e'^, 62 = e-^^ = e^"^ ^ = l,2,...,d. 



whereupon (4^) becomes 

d 

q"(an+e, + fln-ej " (&iq^" + ^2)011 = Aq^ttn, n G Z"'. (4.6) 

The last formula employs standard multi-index notation, e.g. q" = (H^fj^ ' ' ' ■ 
We let formally 



y(t)= J2 anexp|6fq"t| 



and note that 



y'it) = bf anq"exp|6fq"t 



y 



"(t) = hi J2 a„q2"exp 6fq"ti. 



Therefore, multiplying ( [4.6| ) by explftj^q^tl and summing up for n S yields, after 
brief manipulation, the complex functional differential equation 

y"{t) + b2y{t) = 1^ {q^Wiq^'t) + qeviqit)) - Xy'{t)\ . (4.7) 
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Equation (4/7) is of independent interest, being a special case of the equation 



y"{t) + ciy'(t) + C2y{t) = / y(e'^t) dfz{6) 



2n 



where d^u is a complex-valued Borel measure. This, in turn, is similar to the functional 
integro-differential equations of the form 

y"{t) + ciy'it) + C2y{t) = [' y{qt) dr]{q), 



say, where dry is, again, a complex-valued Borel measure. Equations of this kind have 
been considered by the present author, jointly with Yunkang Liu |11], with an emphasis 
on their dynamics and asymptotic behaviour. In the present paper, however, we are 
interested in the spectral problem for ( [4.7] ), and to this end we again expand y in Taylor 
series, y{t) = X)m=o Vrnt^ /fn\. It is easy to affirm by substitution into ( [4.7D the three-term 
recurrence relation 



ym+i = bi^ {2y2(^osaim- X]ym-b2ym~i, m = l,2, (4.8) 




Note a most remarkable phenomenon - although ( [4. 51) is d-dimensional, the index in (|4.8| ) 
lives in Z+! In other words, the dimensionality of the resultant three-term recurrence 
is independent of d - it is, instead, expressed as the number of harmonics in the re- 
currence coefficient. Moreover, inasmuch as (4.8) is more complicated for d > 2 then 
its one-dimensional counterpart (|l.4| ), both recurrences display similar qualitative char- 
acteristics. In particular, we can use the theory of Section |2| to cater for the case of 
ai/vr, a2/TT, . . . , Od/vr being all rational. 

In line with the analysis of Section ^, we let ym (t) = hi yU-ihb2pt), m G Z+, 
whereupon (4.8) becomes 

ym+i{t) = [t + -^ cos aim\yrn{t) -yrn~i{t), m = l,2, (4.9) 



£=1 / 

To recover all solutions of (|4.9| ) we need to consider a linearly independent two-dimensional 
set of solutions. Letting r_i = 0, ro = 1 and sq = 0, si = 1, we recover, similarly to (1.6- 
7), two sequences {rm}meZ+ {sm}meZ+ that span all solutions of (|4.9| ) and such that 
degr^ = "m, degSm = m — 1. In other words, by the Favard theorem both {rm}m.eZ+ 
and {sm+i}mgZ+ ^^'^ OPS and, in line with our analysis of the one-dimensional almost 
Mathieu equation ( |1.1D , we are in position to exploit the theory of orthogonal polynomials. 
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